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PREFACE 

Many  studies  of  one— dimensional  transport  processes 
can  be  found  in  the  literature  of  mathematical  physics. 

However,  the  literature  mostly  ignores  the  complicated 
problem  of  the  determination  of  probability  distributions 
involved  in  neutron  transport  processes.  In  this  Memorandum, 
the  author  begins  a  probabilistic  investigation  of  transport 
processes  and  arrives  at  some  functional  equations  for  generat¬ 
ing  functions  involved  in  the  process.  This  study  evolved 
from  research  in  neutron  transport  theory  sponsored  by  the 
Advanced  Research  Projects  Agency. 
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SUMMARY 

A  probabilistic  discussion  of  one-dimensional  multi¬ 
type  neutron  transport  processes  is  given,  and  the  invariant 
imbedding  technique  is  applied  to  the  generating  functions 
of  the  emitted  neutrons.  A  system  of  simultaneous  partial 
differential  equations  for  these  generating  functions  is 
then  derived.  The  well-known  Ricatti  equation  for  the 
expected  reflected  flux  follows  from  these  differential 
equations. 
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A  PRELIMINARY  INVESTIGATION  OF 
TRANSPORT-BRANCHING  PROCESSES 

1.  INTRODUCTION 

The  purpose  of  this  Memorandum  is  to  begin  a  probabilis¬ 
tic  study  of  neutron  transport  processes,  in  which  neutron 
multiplication  plays  a  role.  We  consider  a  model  of  the 
one-dimensional  neutron  transport  process  that  has  been 
extensively  studied  from  a  nonprobabilistic  viewpoint  by 
Bellman,  Kalaba,  and  Wing  [2,  3],  and  others.  The  essential 
step  of  this  nonprobabilistic  viewpoint  is  to  define  a  vector 
representing  the  expected  value  of  the  neutron  flux  and  to 
regard  this  expectation  vector  as  being  propagated  through 
the  medium.  This  description  of  the  neutron  transport 
process  then  bears  a  strong  resemblance  to  the  physical 
model  of  radiative  transfer,  first  studied  by  Stokes  in 
the  nineteenth  century.  A  brief  history  of  the  study  of 
radiative  transfer  and  neutron  transport  processes,  as  well 
as  an  extensive  list  of  references,  can  be  found  in  Redheffer 
[6] ;  references  concerning  neutron  transport  are  given  in 
Bellman  and  Kalaba  [1].  In  our  probabilistic  investigation 
of  neutron  transport  processes,  we  also  utilize  the  concept 
of  a  branching  process  14].  We  therefore  refer  to  our 
mathematical  model  of  a  one— dimensional  neutron  transport 
process  as  a  tr  an spo r t-br anching  process. 
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In  the  model  of  a  one-dimensional  neutron  transport  pro¬ 
cess  that  we  shall  use  throughout  this  Memorandum,  our  basic 
assumption  is  that  neutrons  are  point  particles  constrained  to 
move  along  a  finite  line  segment.  This  finite  line  segment, 
which  we  shall  refer  to  as  the  rod,  represents  the  inter¬ 
acting  medium  through  which  the  neutrons  travel.  The 
neutrons  travel  through  the  rod,  moving  either  to  the 
right  or  to  the  left,  and  interact  with  the  medium,  occasion¬ 
ally  producing  additional  neutrons  by  nuclear  fissions.  We 
assume  that  the  neutrons  do  not  interact  with  each  other, 
so  that  the  probability  of  a  given  neutron  entering  into 
a  nuclear  fission  or  otherwise  interacting  with  the  medium 
is  independent  of  the  presence  of  other  neutrons.  We 
further  assume  that  these  probabilities  are  independent  of 
time.  In  general,  energy  dependence  may  be  involved,  but 
we  shall  restrict  the  allowable  energies  of  the  neutrons 
to  a  finite  number  of  values,  called  energy  states.  A 
neutron  is  then  completely  described  by  its  position  in 
the  rod,  the  direction  it  is  traveling  ("to  the  left"  or 
"to  the  right"),  and  its  energy  state.  A  neutron  transport 
process  is  initiated  by  neutrons  entering  the  rod  at  the 
left  or  "he  right,  and  when  a  neutron  leaves  the  rod  (at 
the  left  or  the  right),  it  can  no  longer  interact.  The 
models  oi  one-dimensional  neutron  transport  processes 
given  in  [1,  2,  3]  are  similar  to  our  model  and  differ  by 
varying  degrees  of  generalization  and  specialization. 


The 
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description  of  the  physical  parameters  needed  to  describe 
the  process  is  given  in  Sec.  5  of  this  Memorandum. 

We  are  interested  in  the  problem  of  determining  the 
probability  distributions  of  the  neutrons  leaving  the  rod 
in  neutron  transport  processes  initiated  by  given  beams 
of  neutrons.  One  could  consider  the  neutron  transport 
process  as  a  continuous-time  stochastic  processes,  but 
such  an  approach  does  not  yield  explicit  solutions  to 
the  particular  problem  that  we  are  considering — determining 
the  probabilities  that  given  numbers  of  neutrons  leave  the 
medium  (rod)  throughout  the  process.  The  general  approach 
that  we  use  is  to  consider  the  probability  distributions 
of  the  neutrons  leaving  the  medium  as  functions  of  the  ini¬ 
tial  and  end  points  of  the  rod  (which  we  consider  as  imbedded 
in  the  real  line).  This  technique,  for  which  Bellman  coined 
the  term  invariant  imbedding  [1,  2,  3],  has  been  used  by 
Bellman  and  others  in  the  study  of  neutron  transport  proces¬ 
ses  as  well  as  in  a  variety  of  other  problems.  We  first 
consider  the  problem  of  determining  joint  probability  dis¬ 
tributions  of  the  neutrons  "reflected"  and  "transmitted" 
through  a  rod  composed  of  two  adjoining  rods,  each  with 
known  probabilities  for  reflection  and  transmission.  The 
approach  that  we  use  is  similar  to  the  method  employed  by 
Mycielski  and  Paszkowski  [5]  in  discussing  a  simple  particle 
transport  process,  which  is  a  special  case  of  the  process 
we  consider.  A  direct  application  of  the  Myci el ski-Pa szkowski 
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method  to  our  problem  would  result  in  immense  combinatorial 
difficulties  unless  a  proper  choice  of  random  variables  is 
made.  We  use  a  set  of  random  variables  that  describe  a 
multitype  (Galton-Watson)  branching  process  (4],  and  we 
are  thus  able  to  obtain  simple  functional  equations  for 
the  generating  functions  describing  the  reflected  and  trans¬ 
mitted  neutrons.  Finally,  by  using  a  limiting  process,  we 
obtain  partial  differential  equations  for  these  generating 
functions. 
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2 .  BRANCHING  PROCESSES 

We  begin  by  summarizing  the  terminology  and  notation 
of  branching  processes  [4]  that  we  use  throughout  this 
Memorandum.  A  branching  process  consists  of  a  collection 
of  objects — e.g.,  neutrons  in  a  chain  reaction — each  of 
which  can  produce  similar  additional  objects.  If  we  con¬ 
sider  an  initial  collection  of  objects  as  belonging  to 
"generation  0,"  then  the  collection  of  all  objects  produced 
by  these  objects  are  considered  as  objects  of  "generation  1." 
In  general,  the  objects  of  generation  n  produce  objects  of 
generation  n  +  1.  Since  the  probability  distribution  of 
the  objects  of  a  given  generation  depends  only  on  :he  dis¬ 
tribution  of  objects  of  the  immediately  preceding  generation, 
we  thus  have  a  Markov  chain 


where  Zn  describes  the  objects  of  the  n-th  generation,  for 
nonnegative  n.  Furthermore,  the  distribution  of  the  descen¬ 
dants  of  a  given  object  of  generation  n  is  independent  of 
the  distribution  of  decendants  of  the  other  objects  of 
generation  n. 

We  shall  consider  branching  processes  involving  a  finite 
number  of  types  of  objects,  and  we  number  the  types  from  1 
to  m.  The  state  of  the  process  in  generation  n,  represented 
by  Zn  in  the  Markov  chain,  is  described  by  an  m-dimensional 
vector  with  nonnegative  integer  components,  (z^,  ...,  z^)  , 
where  z^  is  the  number  of  objects  of  type  i  in  generation  n. 
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We  can  thus  consider  Zn  as  a  vector-valued  random  variable, 
which  we  shall  refer  to  as  a  random  vector.  The  generating 
function  g^  of  the  random  vector  Zn  is  defined  by 


(2.1) 


where 


p(Zn  -  (k,,  k 


2* 


k  )) 
m 


for  n  -  0,  1,  2,  -  We  consider  the  domain  of  g,  x  to 

m  (n) 

be  D  ,  where  D  is  the  closed  unit  disk  in  the  complex  plane, 

and  we  note  that  g(n)  :  Dm  D  and  that  g(n)(l,  ...,  1)  -  1. 

We  quote  now  the  basic  result  of  Kolmogorov  and  Dmitriev 

(see  Harris  [4],  p.  36)  for  branching  processes.  Let 

ak  ...k  (f°r  i  "  I?  •••,  m;  k,  ,  ...,  k  *0,  1,  2,  ...) 

represent  the  probability  that  a  particle  of  type  i  will 

produce  exactly  k^  particles  of  type  j  for  j  -  1,  .  .  .  ,  m. 

We  define  the  function  f :  Dm  -*  Dm  by 

f(s)  -  (f^s),  ....  £nl(s))  , 


where 


fl(s) 


} 


(2.2) 


~  “  W’  •  •  •  '  sm^  • 
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The  resulting  relation  describing  the  generating  functions 
fg(n)^  is  given  by 

(2’3)  S(n)(s)  -  8(0)(f(n)(£))  • 


where  f^n^  is  the  n— th  iterate  of  £  (defined  recursively  by 

~(n+l)  ’  £(0)(£)  "  $)  ■  We  cali  £  the  branching 

generating  function  of  the  branching  process  fZn1, 

A./ 

Some  notation  that  we  find  convenient  to  use  throughout 

this  Memorandum  is  summarized  below.  We  denote  vector 

quantities  by  an  underscored  and  write  s  -  (s, ,  .  ..,  s  ), 
Ill  ~  1  m 

£  *  (s^  •  •.,  sm),  etc.  The  symbol  1.  represents  the 

vector  (1,  .  ..,  1).  We  also  use  this  convention  for  sub¬ 
scripts,  and  the  symbol  ?  means  the  sum  over  all  nonnegative 
values  of  jr  jr  ....  and  jm  (where  £  =  (jr  . jj). 

For  example,  Eq.  (2.1)  becomes 


g 


(n) 


We  also  utilize  the  concept  of  the  conditional  generat¬ 
ing  function  of  a  vector-valued  random  variable.  If  W  is 
an  m-dimensional  random  vector,  and  (X,  is  an  event,  then 
the  conditional  generating  function  c o  of  (W| &  )  is  defined  by 


k, 

®(s)  -  2  P(W  -  WCL  )  s, 1. . .  s  m 
^  ~  ~  i  m 


(2.4) 
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3.  THE  TRANSPORT  GENERATING  FUNCTIONS 

We  now  consider  the  one-dimensional  neutron  transport 
process  described  in  Sec.  1.  Let  the  process  take  place 
on  a  rod  extending  from  0  to  x  on  the  real  line  and  assume 
that  there  are  m  distinct  energy  states  numbered  from  1  to 
m.  In  order  to  allow  a  simple  mathematical  solution  of 
this  one-dimensional  transport-branching  process,  we  shall 
consider  the  following  variables: 


-  the  number  of  neutrons  of  energy  state  i  leaving 
the  reac.or  at  the  left, 

Y,  ■  the  number  of  neutrons  of  energy  state  i  leaving 
the  reactor  at  the  right, 

i  *  1,  . .  .  ,  m. 


We  also  define  the  quantities: 


*  the  number  of  neutrons  of  energy  state  i  incident 
at  the  left  (traveling  toward  the  right), 

Vi  -  the  number  of  neutrons  of  energy  state  i  incident 
at  the  right  (traveling  toward  the  left), 

i  *  1,  »  »  •  ,  m. 


(See  Fig.  1. ) 


Fig.  1 
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We  then  define  the  random  vectors 


(X, ,  . . . ,  xm) , 


I  -  (Yr  ....  Ym), 
U  -  •••>  V-  l-  <vr  •••>  V  ), 


and 


5  "  (S’  X)  ■  (*1 . V  Yr  ....  Ym)  . 

We  also  let  E1  represent  the  "i-th  unit  vector"  with  the 
1  th  coordinate  equal  to  one,  and  other  coordinates  zero. 

Since  the  number  of  neutrons  emitted  from  the  reactor 
depends  on  the  incident  neutron  beams,  the  generating 
function  for  Z  depends  on  the  assumptions  about  the  values 
of  U  and  V.  It  is  useful  to  define 


Y 


£  -  the  conditional  generating  function  of 
( Z |  U  «  0,  V  =  E1)  , 


"  the  conditional  generating  function  of 


(Z|U  -  E1,  V  -  0)  , 


for  i  -  1,  . .  .  , 


m. 


We  also  define  the  vector  generating  functions  Y^  and  Y^  by 
T<-(s)  -  (if-(s) . ^(s)),  1-1,2. 

Assume  that  U  and  V  are  random  vectors,  and  let 
2m 

g:  D  -  D  be  the  generating  function  of  the  2m-dimensional 
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vector  (V^  •••>  Vm,  U^,  .  ..,  Um).  Since  the  neutrons  do 
not  interact  with  each  other,  it  is  easily  seen  that  the 
generating  function  for  Z  is  given  by 

gQ^s),  r2(s))  . 

I*SJ  <»v 

Our  eventual  goal  is  to  determine  the  generating 
1  2 

functions  and  ,  which  we  call  the  transport  generating 
functions .  Following  the  Mycielski-Paszkowski  method  [5] 
of  considering  multiple  "reflections"  and  "transmissions," 
we  first  determine  the  transport  generating  functions  for 
a  rod  composed  of  two  segments  in  terms  of  the  transport 
generating  functions  of  the  segments. 

Consider  the  composite  rod  a3  made  up  of  the  two  rods 
a  and  3  as  in  Fig.  2.  We  let  Y^,  and  (/  -  1,  2) 
be  the  transport  generating  functions  for  the  neutron 
transport  processes  on  the  rod  a,  the  rod  3,  and  the  rod  a3, 
respectively.  Let  i  be  a  fixed  integer  between  1  and  m, 
and  consider  the  transport-branching  process  initiated  by 
one  neutron  of  energy  state  i  incident  upon  the  rod  a3  at 
the  right  (and  no  neutrons  incident  at  the  left).  Let  X 
and  Y  be  the  random  vectors  representing  the  neutrons 
emerging  from  the  rod  oc3,  as  defined  above.  The  generating 
function  for  Z  *»  (X,  Y)  is  thus  given  by  the  i— th  component 
of  •  This  process  is  to  be  viewed  as  consisting  of 
multiple  "reflections"  and  "transmissions"  through  the 
segments  a  and  3  (see  Fig.  2).  We  let  (V°,  R°)  represent 
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the  neutrons  emerging  directly  from  P  with  the  single 
neutron  of  energy  state  i  incident  at  the  right.  Then  we 
let  (T^,  U1)  represent  the  neutrons  emerging  directly 
from  a  with  incident  neutron  beam  represented  by  V°  at 
the  right.  In  general,  (Vn,  Rn)  represents  the  neutrons 
emerging  directly  from  p  with  incident  neutron  beam  11°  1 
at  the  left;  (Tn"^,  Un"*"^)  represents  the  neutrons  emerging 
directly  from  a  with  incident  neutron  beam  Vn  at  the  right. 
(Vn  and  Rn  are  considered  as  0  for  n  odd;  Tn  and  Un  are 
considered  as  0  for  n  even.)  We  define 

Xn  *  Z  Tk  , 

~  k«0~ 

(3.1) 


,n 


n  l 
•  2  R 

k-0~ 


This  complicated  description  can  be  greatly  simplified 
by  considering  the  branching  process  involving  4m  objects, 
where  the  object  of  type  i  represents 

a  neutron  of  the  energy  state  i  emerging  from 
ap  to  the  left,  for  1  <  i  <  m; 
a  neutron  of  energy  state  i  -  m  emerging  from 
J  ap  to  the  right,  for  m  +1  <  i  <  2m; 

a  neutron  of  energy  state  i  -  2m  emerging  from 
P  to  the  left,  for  2m  +1  <  i  <  3m; 
a  neutron  of  energy  state  i  -  3m  emerging  from 
a  to  the  right,  for  3m  +  1  <  i  <  4m. 
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This  branching  process  is  described  by  the  Markov 

chain 

n°  n1  nn 

•••*;£»  •••  » 

where  the  4m -dimensional  random  vector 


(3.2) 


S"  -  (X".  Vn,  Un)  ,  for  n  -  0,  l,  2,  ...  . 


Let  be  the  generating  function  of  gn  and  let  f :  D4”1  -  D4™ 

bethe  branching  generating  function  for  the  branching  process 

{Qn}. 

Writing 

(3.3)  f  -  (f1,  f2,  f3,  f4)  , 

^  ~ 

where  f*':  D4™  -•  Dm,  /  *  1,  . ..,  4,  the  reader  can  easily 
see  from  the  above  description  that 

£l(sl >  ?2>  s3,  _s4)  =  s1  , 


(3.4) 


f2(s1,  s2,  s3,  s4)  «  s2  , 
f3(s\  s2,  s3,  s4)  -  As1,  s4)  , 
f4(sx,  s2,  s3,  s4)  «Y?(s3,  s2)  . 

Writing  f(n)  -  (£(n)>  £(n)>  |(n) »  £(n)^  we  have  the  recursive 
formulas : 


1 

I 
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(3.5) 


'(n) 
2 

(n) 

3 

(n) 

~(n) 


li^1'  ^(n-2)’  S2))  • 

l|(£(n-l)-  l2)  £(n-2))’  £2) 


The  Markov  chain  Qn  will  converge,  providing  Vn 
approaches  zero  (or  equivalently,  Un  -*()).  Consider  the 
Markov  chain 


9 


which  also  describes  a  branching  process  with  branching 
generating  function  given  by 

(3.6)  f  (s)  -Y^U,  Y?(s,  1))  . 


If  we  assume  that  the  Jacobian  of  f  (evaluated  at  1)  is 

~v  - — ' 

small  enough,  then  Vn  -  0  with  probability  one;  and  it 
follows  that  Un  -*  0  with  probability  one,  and  2°  converges 
to  a  finite  value  with  probability  one. 

We  henceforth  assume  the  convergence  (with  probability 
one)  of  Qn.  This  means  that  the  probability  of  a  "chain 
reaction"  (producing  an  infinite  number  of  neutrons)  occur¬ 
ring  is  zero.  Then  the  generating  functions  and 

8(n)  *  S(0)£(n)  conver8e>  and  lim  g(n)  is  the  generating 
function  of  lim  Qn.  From  the  definition  of  Qn,  it  is  clear 
that 


d 
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(3. 7)  lim  Qn  -  (x,  Y,  0,  0) 
n-»® 


Let 


(3.8) 


fCn)^1.  £2-  0,  0)  =  ^(s1,  s2)  , 

n-»  ~(n)^~lj  £2,  2’  2>  -J)V-  £2)  • 


Let  cp^  and  cp^,  be  the  i-th  components  of  Y^  and  y^d, 


a3 


respectively.  Then 


(o)  ■  Vs3>  £2)  > 


and  therefore, 


S(n)  =  ®g^£(n)  ’  £2)  • 

The  generating  function  for  lim  gn  -  (X,  Y,  0,  0)  is  given  by 
8  "  lim  8(n)  “  s2),  s2)  . 


Therefore, 


^ap^1,  Z2)  “  S^1,  s2,  0,  0)  «  CDpC^Cs1,  s2),  s2)  , 


/V  ^ 


from  which  we  conclude  that 


(3.9) 


lip  ■  ZpCu1.  s2)  . 


'  11  • 
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By  identical  reasoning,  or  by  appealing  to  symmetry, 
we  obtain  our  second  relation 

(3-10>  1^(3  "  Ia(£lj  ~2)  ' 

1  2 

The  functions  tt  and  tt  ,  which  are  the  limits  of  the  sequences 
given  by  Eqs.  (3.8),  satisfy  the  implicit  equations, 


(3.11) 


TT1  =  Y^Cs1,  tt2)  =  Y^s1,  v 2 ( tt1  ,  S2)) 


n2  s2)  =  Xf^Y^s1,  n2),  s2) 
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_ DIFFERENTIAL.EQUATIQNS  OF  INVARIANT  IMBEDDING 

We  now  use  the  invariant  imbedding  approach  combined 

with  the  functional  equations  of  Sec.  3  to  obtain  systems 

of  partial  differential  equations  describing  the  transport 

1  o 

generating  functions  Y  and  Y  .  One  should  note  that  these 
partial  differential  equations  could  also  be  obtained,  using 
the  techniques  of  this  Memorandum,  without  recourse  to  the 
results  of  Sec.  3. 

Let  us  first  consider  a  neutron  transport  process 

taking  place  on  a  given  rod  [0,  L] .  To  utilize  the  invariant 

imbedding  method,  we  consider  the  family  of  segments  [a,  b] 

of  the  rod  [0,  L]  ,  and  we  let  yJ  ,  ,  and  v?  be  the 

~ia,Dj  ~la,bj 

transport  generating  functions  defined  in  Sec.  3  for  the 
transport  process  taking  place  on  the  (isolated)  rod 
segment  [a,  b]  (for  0  <  a  <  b  <  L) . 

Write 


lfo,x](£>  "  ~4(x’  s)’  1  ■  1’  2  ’ 

where 

3  -  (s1,  s2)  e  D2™  , 

and  consider  the  family  of  generating  functions  {Y^(x,  s)} 
parametrized  by  the  ral  variable  x.  In  order  to  obtain 
differential  equations  for  Y^x,  s),  we  consider  the  rod 
[0,  x  +  h]  as  composed  of  two  segments  [0,  x]  and  [x,  x  +  h] . 
(See  Fig.  3 . ) 
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a 

a 


0 


->r 


0 


x  x  +  h 


Fig.  3 


Letting 


<4-1)  3K  ifx.x+h]^ 


h=0 


=  u/(x,  s),  l  ~  l,  2  , 


we  have 


(4.2) 


Y[XjX+h](?1<  i>2)  "  s/  +  U)'(x,  s)h  +  o(h)  ,  <,  =  1,  2 


1  2 

(We  shall  assume  that  uj  and  are  continuously  differentiable.) 

Applying  the  results  of  Sec.  3  with  a  =  [0,  x]  and 
p  =  [x,  x  +  h]  (see  Fig.  3),  we  have 


Ax  +  h,  s1,  ,2)  =  ifx.x+hjAC?1’  ?2).  s2) 


(4.3) 


=  TT^(sS  s2)  +  u)^(x,  tY(sY  s2),  s2)h+o(h) 


Y2(x  +  h,  s1,  s2)  «  Y2(x,  s1,  tt2(s1,  s2))  . 


Applying  Eqs.  (3.8),  we  find  that  the  functions  rr  and 

2 

n  in  the  above  expression  are  given  by 
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ar 


(4.4) 


~  ~  ^  “  1  £  5  h  ?  ^  — K  +  o,(h), 

j  »sf 


.2,1  .2 


'£'  +  JS^Cx.  I1,  s2)h  +  o(h)  . 


We  readily  obtain  from  Eqs.  (4.3)  and  (4.4)  the  basic  system 

o_f  partial  differential  equations  for  the  transport  generating 
functions , 


k  il<**  i1-  52) 


*Y 


t  J:1-  ^2)  71  +  1 1<  £2)  . 


(4.5) 


lie  I2(x>  s1,  s2) 


^Sj 


BY' 


2*  w.(x,  _Y  ,  ,s2)  — .  . 
J  2  Zsj 


The  initial  conditions  are  given  by 


T1(0,  s1,  s2)  -  s1 


(4.6) 


12(0,  sl,  s2)  «=  s2 


Note  that  the  first  equation  of  (4.5)  does  not  involve  Y2, 
so  it  can  be  solved  independently  to  give  the  value  of  Y1. 

We  can  also  regard  these  partial  differential  equations  as 
an  nfinite  system  of  ordinary  differential  equations  for 

the  probabilities  of  various  combinations  of  particles  being 
reflected  and  transmitted. 
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If  we  are  interested  only  in  the  "reflection"  probabili¬ 
ties,  these  basic  equations  can  be  simplified  by  letting 

(4.7)  $(x,  s)  -  Y^Cx,  1,  s)  . 

^  ^ 

The  reader  should  observe  that  $.(x,  s)  is  the  conditional 
generating  function  of  (Yl  U  *  0,  V  «=  E1)  in  the  transport 
process  taking  place  on  the  rod  [0,  x] .  The  first  equation 
of  (4.5)  becomes 

%  £(x> 

(4.8) 

*(o,  s) 

This  last  result  is  analogous  to  the  matrix  Ricatti  equation 
for  the  expected  reflected  neutron  flux  [1].  In  fact,  the 
Ricatti  equation  for  this  reflection  function  follows  from 
Eq.  (4-8),  as  is  shown  at  the  end  of  the  next  section. 


m 

Z 


2 

ju  .  (  x , 
J 


s)  p—  $(x,  s)  + 

~  j  ^  ^ 


If 

a-  (x. 


$,  s) 
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5.  EQUATIONS  FOR  NEUTRON  TRANSPORT 

The  differential  equations  of  Sec.  4  are  dependent  on 
Eq.  (4.2)  (or,  equivalently,  on  Eq.  (4.1)),  which  defines 

I  r\ 

the  key  parameters,  a;  and  tu  ,  of  an  abstract  transport¬ 
branching  process.  We  shall  now  determine  uu^  and  in 
terms  of  some  physical  parameters  of  a  neutron  transport 
process.  Let 


^(x)  =  The  reciprocal  of  the  mean  free  path  of  a 

neutron  of  energy  state  i  at  the  point  x  in 
the  rod, 

Pr°bability  that  a  nuclear  fission  initiated 
by  a  neutron  of  state  i  at  the  point  x  will 
produce  exactly  j  neutrons  of  state  u  travel¬ 
ing  to  the  left  and  k^  neutrons  of  state  u 
traveling  to  the  right,  for  u  -  1,  . . .  ,  m. 


(Absorption  of  neutrons  can  be  taken  into  account  by  con¬ 
sidering  absorption  as  a  "nuclear  fission"  in  which  no 
neutrons  are  produced. )  We  then  define  the  generating 
functions 


(5.1)  Q^(x,  s,  t) 

1  cv 


(x)s 


for  i  ■  1, 


m. 


Consider  a  neutron  fission  process  in  the  rod  [x,  x  +  h] . 
A  neutron  of  energy  state  i  entering  the  rod  has  the  prob¬ 
ability  fTi(x)h  +  o(h)  of  causing  a  nuclear  fission,  and  the 
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probability  that  this  neutron  starts  a  chain  reaction  result¬ 
ing  in  two  or  more  fissions  is  o(h).  Letting  represent 
the  i-th  component  of  y|x  x+h]0  =  1>  2;  i  =  1,  .  .  .  ,  m)  ,  we 
therefore  conclude  that 


[1  -  rT^(x)h  +  o(h)]s£  (no  fission), 
+  [a^(x)h  +  o(h)]Q^(x,  s^,  £2)  (one  fission), 
+  o(h)  (two  or  more  fissions), 


or,  more  simply, 


(5.2)  cd^s1,  s2)  =  s*  +  a- (x)[Q.(x,  s1,  s2)  -  sf]h  +  o(h) 


Comparing  Eq.  (5.2)  with  Eq.  (4.2),  we  obtain  our 
desired  result: 

(5.3)  w'(x,  s1,  s2)  =  a.(x)[Q.(x,  s1,  s2)  -  sU  , 

1  (V;  ^  X  1  «n/  X 

for  i  =  1,  2;  i  =  1,  .  . .  ,  m. 
With  only  one  energy  state  (m  =  1),  Eqs.  (4.5)  become 

1  BY1 

^CxT  ^x 

1  ^Y2 

oTx)  "5x 


=  [Q(x,  Y1,  s2)  -  s2]  -Y-^  +  Q(x,  Y1,  s2)  -  Y1 


[Q(x,  Y1,  s2)  -  s2]  — 


(5.4) 
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The  neutron  fission  process  in  which  each  fission 
produces  two  neutrons,  one  traveling  to  the  left  and  one 
traveling  to  the  right,  is  often  discussed  in  the  literature 
[1,2,3].  In  th  is  case, 

(5.5)  Q(x,  s^,  s2)  =  sV  ; 


the  transport  generating  functions  are  given  by 

oCxT  -  s2)  ^  +  s V  -  y1  , 


(5.6) 


•Y 


?u,l  „  2\  ftY‘ 


"(s^-s^) 


?)s 


’f1(0,  S1,  s2)  =  s1,  y2(0,  S1,  s2)  =  s2  , 


and  the  generating  function  for  the  reflected  neutrons  is 
given  by 


1  -H  /_.  N 
^TxT  "  (s$  ~  s)  B?  +  s*  “  «  > 

(5.7) 

*(0*  s)«  1  . 

This  special  case  is  discussed  by  Bellman,  Kalaba, 

and  Wing  [2],  who  arrived  at  Eq.  (5.7)  using  elementary 
methods. 

Another  simple  example  of  a  transport-branching  process 
is  the  continuous  one-dimensional  random  walk.  In  a  one¬ 
dimensional  random  walk,  a  particle  enters  the  rod  and 


r 
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random  ly  reverses  its  direction  during  its  journey  until 
it  leaves  the  rod  either  to  the  left  or  to  the  right.  Let 


o/x) 

"r(x) 


the  reciprocal  of  the  mean— free  path  of  the 
particle  when  it  is  traveling  to  the 
fief  t  ^ 
right) 


at  the  point  x. 


Thus  the  particle  has  the  probability  ar(x)h  +  o(h)  of 
reversing  its  direction  while  traveling  to  the  right  from 
x  to  x  +  h,  and  similarly  while  traveling  to  the  left. 

The  reader  can  then  readily  see  that 

U)1(x,  s,  t)  =  a  (x)(t  -  s)  , 

(5.8) 

a'  (Xj  s ,  t)  =  ar(x)(s  -  t)  . 


Equation  (4.8)  then  becomes 

||  *=  "r(>0($  -  s)  ||  +  at(x)(s  -  $)  , 

(5.9) 

$(0,  s)  *  1  . 

The  solution  to  (5.9)  is  given  by 

$(x,  s)  =  1  -  p(x)  +  p(x)s  , 

where  p(x)  is  the  solution  of  the  Ricatti  equation 

2 

p'  =  c,  ~  (r^  +  ar)p  +  arp  , 


(5.10) 


p(0)  =  0  . 
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The  quantity  p(x)  is  the  probability  that  a  particle  enter¬ 
ing  the  rod  [0,  xj  at  x  will  leave  the  rod  at  x,  and  this 
probability  is  given  by  the  well-known  equation  (5.10). 

The  well-known  matrix  differential  equations  [1,  6]  of 
a  neutron  transport  process  are  easily  derived  from  our 
general  results.  Let  R^(x)  represent  the  expected  number 
of  neutrons  of  state  i  leaving  the  rod  [0,  x]  at  the  right, 
given  that  one  neutron  of  state  j  enters  the  rod  at  the 
right.  Then 


(5.11) 


If  the  probability  that  a  finite  number  of  neutrons 
are  produced  is  equal  to  one  (i.e.,  the  probability  of  the 
occurrence  of  a  "chain  reaction"  is  zero),  it  follows  that 


$t(x,  1)  =  1  , 

w£(x,  1)  »  0  ,  i  -  1,  2  , 


i 


1, 


•  •  y  m  • 


Differentiation  of  Eq.  (4.8)  utilizing  the  above 
identities  yields  the  familiar  matrix  Ricatti  equation  [1] 
for  the  expected  reflected  neutron  flux. 
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(5.12) 


d  R 

<35  Rij 


+  2  R 


ik 


+ 


2 

k,l 


9 


1  2 

where  the  partial  derivatives  of  uu  and  uu  are  evaluated 
1  2 

at  s  ■  s  =1.  The  Jacobian  matrices 

<v  ~ 


/V  Arf  AW 


1  2 

(evaluated  at  s  =  s  =  1.)  appearing  in  Eq.  (5.12)  have  the 
following  physical  significance.  Consider  a  neutron  trans¬ 
port  process  taking  place  on  the  isolated  rod  segment 
P  «  [x,  x  +  h] ,  in  which  one  neutron  of  energy  state  j 
enters  at  the  right  (and  no  neutrons  enter  at  the  left). 
Then  the  expected  number  of  neutrons  of  energy  state  i 
leaving  p  at  the  right  is  given  by 


dllU 

—4  (x,  1,  1)  +  o(h)  , 


and  the  expected  number  of  neutrons  of  state  i  leaving  at 
the  left  is  given  by 


f 

ij  “  t~t  b  V  +  °(h)  • 

J  OS^ 
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2  2 

B(JU  Buu 

The  matrices  — j-  and  — *  are  similarly  interpreted. 
Bs  Bs^ 


vVe  can  similar ly  define  T^j(x)  to  be  the  expected 
number  of  neutrons  of  state  i  leaving  the  rod  [0,  x]  at 

the  left,  given  that  one  neutron  of  state  j  enters  the  rod 
at  the  right.  Then 

(5.13)  ^.(x)  “  ~4  (x>  1)  , 

and  again  assuming  that  the  probability  of  a  "chain  reaction" 
is  zero,  we  similarly  obtain  the  differential  equation  for 
the  transmission  matrix 


(5.14) 


d  T 

as  Tu 


OUU, 

s  Tu<rt +  ?  -T  <W 


k  .s1 


Finally,  we  remark  that  the  algebraic  equations  for 
nonprobabilistic  transport  theory  found  in  Redheffer  [6] 
are  also  easily  obtained  from  the  results  of  Sec.  3  by 
taking  expectations. 
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